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Abstract 

We calculate the cross section of charm photoproduction at HERA collider in the frame- 
work of the fcy-factorization QCD approach. Our analysis cover the inclusive charm pro- 
duction as well as charm and associated jet production processes. Both photon-gluon and 
gluon-gluon fusion mechanisms are taken into account. The unintegrated gluon densities 
in a proton and in a photon obtained from the full CCFM, from unified BFKL-DGLAP 
evolution equations as well as from the Kimber-Martin-Ryskin prescription are used. Our 
theoretical results are compared with the recent experimental data taken by the HI and 
ZEUS collaborations at HERA. Special attention is put on the specific angular correlations 
which can provide unique information about non-collinear gluon evolution dynamics. 



1 Introduction 

The charmed quark production in electron-proton collisions at HERA is a subject of 
intensive study from both theoretical and experimental points of view [1-6]. The value of 
charm mass m c provides a hard scale which allows perturbative QCD (pQCD) to be applied. 
The production dynamics is governed by the photon-gluon 75 — > cc or gluon-gluon fusion 
gg — > cc (direct and resolved photon contributions, respectively) and therefore cross sections 
of such processes are sensitive to the gluon content of a proton and of a photon. Very 
recently, the HI and ZEUS collaborations have presented important experimental data [5, 
6] on the charm photoproduction at HERA. In [6] the data sample five times larger than 
in previous analysis [1] has been used. Differential cross sections are determined for events 
with a D*^ meson (inclusive D*^ production) and for events with a D*^ meson and one 
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or two hadronic jets. In the ZEUS analysis [5] the differential inclusive jet cross sections 
for events containing a D* ± meson have been measured and specific angular correlations in 
the and dijet associated photoproduction have been studied. A comparison of these 
measurements with the next-to-leading (NLO) pQCD calculations shows [5, 6] that NLO 
pQCD has some marked problems in description of experimental data. In particular, the 
significant differences are observed [5] in the shape and in the normalization of most of the 
distributions between data and theoretical predictions. On the other hand recently much 
progress has been made towards a global understanding of the /^-factorization [7, 8] (or 
semihard [9, 10]) approach by working out this picture for several heavy quark and prompt 
photon production processes at HERA and Tevatron [11-15]. 

The /^-factorization approach is based on the familiar Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) [16] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [17] gluon evolution. In this 
way, the large logarithmic terms proportional to In 1/x are summed up to all orders of 
perturbation theory (in the leading logarithmic approximation). It is in contrast with the 
popular Dokshitzer-Gribov-Lipatov-Altarelli-Parizi (DGLAP) [18] strategy where only large 
logarithmic terms proportional to ln/i 2 are taken into account. The basic dynamical quan- 
tity of the /cy-factorization approach is the so-called unintegrated (i.e. k^-dependent) gluon 
distribution A(x, kf,, /i 2 ) which determines the probability to find a gluon carrying the lon- 
gitudinal momentum fraction x and the transverse momentum k T at the probing scale fj 2 . 
The unintegrated gluon distribution can be obtained from the analytical or numerical so- 
lution of the BFKL or CCFM evolution equations. Similar to DGLAP, to calculate the 
cross sections of any physical process the unintegrated gluon density A(x, k 2 ^, /i 2 ) has to be 
convoluted [7-10] with the relevant partonic cross section. But as the virtualities of the 
propagating gluons are no longer ordered, the partonic cross section has to be taken off 
mass shell (kj-dependent). It is in clear contrast with the usual DGLAP scheme (so-called 
collinear factorization). Since gluons in the initial state are not on-shell and are charac- 
terized by virtual masses (proportional to their transverse momentum), it also assumes a 
modification of their polarization density matrix [7, 8]. In particular, the polarization vector 
of a gluon is no longer purely transversal, but acquires an admixture of longitudinal and 
time-like components. Other important properties of the //^-factorization formalism are the 
additional contribution to the cross sections due to the integration over the k 2 - region above 
/i 2 and the broadening of the transverse momentum distributions due to extra transverse 
momentum of the colliding partons. 

Concerning the theoretical treatment of charm production in the framework of standard 
(collinear) QCD, two types of NLO calculations are available for comparison with the recent 
HI and ZEUS experimental data [5, 6]. The traditional massive charm approach [19] or 
fixed-flavour-number scheme (FFNS) assumes that light quarks are the only active flavours 
in the structure functions of the proton and photon, so that charmed quarks are produced 
only in the hard process. This scheme should be reliable when the transverse momentum pt 
of the charmed quarks is of similar size compared to m c and breaks down for ^> m c . It 
is because of presence of collinear singularities which having the form a s ln(p^/m 2 ). In the 
massless or zero-mass variable-flavour-number scheme (ZMVFNS) [20] charmed quarks are 
treated as an additional active flavours (massless partons). This approach is applicable at 
high transverse momenta pr m c . To completeness, we should also mention the general- 
mass variable-flavour-number scheme (GMVFNS) [21] which combines the massless and 
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the massive scheme. Note that the massless charm calculations take into account charm 
excitation processes and thus predict a larger resolved component in comparison with the 
massive calculations. However, both massless and massive approaches underestimate [6] the 
measured cross section of the inclusive photoproduction in the intermediate transverse 
momentum p T and forward pseudo-rapidity i] regions. An agreement between the theoretical 
and experimental results can only be achieved using some extreme parameter values. In 
particular, in the NLO massive scheme a very low charmed quark mass m c = 1.2 GeV was 
required [19]. But even within this set of parameters, the shapes of the transverse 
momentum and pseudo-rapidity distributions cannot be said well reproduced. Some better 
agreement between the massless scheme and the measured px (though not if) spectrum was 
achieved using special assumptions on the c — > D* fragmentation. The similar situation is 
observed in the case of semi-inclusive charm production: the data tends to agree with the 
upper bound of the NLO calculation [5]. However, the difference between the results of 
NLO calculation observed in shape of inclusive differential cross section as function of the 
pseudo-rapidity i] D *± is not seen in the semi-inclusive cross section as a function of rf ot , and 
the shape of the data is well described by the NLO QCD predictions. At the same time 
the experimentally obtained dijet angular correlations [5] show a large deviation from the 
massive NLO QCD predictions, specially for the resolved-enriched sample. In general, it 
was concluded [6] that for the precise description of the charm photoproduction higher-order 
corrections or implementation of additional parton showers in current NLO calculations are 
needed. 

In our previous paper [11] the ability of the /c T -factorization approach to reproduce the 
recent experimental data for the D* ± and dijet associated photoproduction (including the 
specific angular correlations between the hadronic jets in final state) taken by the ZEUS 
collaboration has been investigated. It was demonstrated [11, 22, 23] that the leading-order 
off-shell matrix elements of the photon-gluon fusion \M\ 2 { r )g* — > cc) combined with the non- 
collinear evolution of gluon densities in a proton A(x, k 2 -, /z 2 ) effectively simulate contribution 
from the charmed quark excitation subprocess eg — > eg. Next, some c — c correlations in 
high energy 7p scattering have been studied [24] and comparisons to the recently measured 
data of the FOCUS collaboration at Fermilab were made. In particular, it was shown that 
the analysis of the kinematical correlations of charmed quarks opens new possibilities for 
verifying models of non-collinear gluon evolution. 

In the present paper we will analyse the recent HI and ZEUS data [5, 6] using the 
£; T -factorization approach of QCD. Mostly we will concentrate on the D*^ and single jet 
associated photoproduction. It is because these processes have not been studied yet in the 
framework of /^-factorization and the comparison to the experimental data [6] were made 
in the framework of MC generator Cascade [25] only. We investigate the different produc- 
tion rates (calculated in a number of different kinematical regions) and make a systematic 
comparison of our predictions to the recent HI and ZEUS data. Special attention will be 
drawn to the specific angular correlations in associated D* ± and jet production since these 
correlations are sensitive to the transverse momentum of the partons incoming to the hard 
scattering process and therefore sensitive to the details of the non-collinear gluon evolution. 
Additional motivation of this study is the fact that /c^-factorization approach automatically 
incorporates the main part of the standard (collinear) high-order corrections [7-10]. Our 
consideration will be based on the leading-order off-shell matrix elements of the photon- 
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gluon and gluon-gluon fusion subprocesses (the direct and resolved photon contributions, 
respectively) which have been calculated in our previous papers [12, 26]. In the numerical 
calculations we will test the different sets of unintegrated gluon distributions in a proton and 
in a photon which are obtained from the full CCFM [27], from the unified BFKL-DGLAP 
evolution equations [28] and from the conventional (DGLAP-based) quark and gluon densi- 
ties. In the last case we will use the so-called Kimber-Martin-Ryskin (KMR) [29] approach. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the /cy-factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Finally, in Section 4, we 
give some conclusions. 



2 Theoretical framework 



2.1 Kinematics 



We start from the gluon-gluon fusion subprocess. Let p e and p p be the four-momenta of 
the initial electron and proton, k\ and k 2 the four-momenta of the incoming off-shell gluons, 
and p c and p 5 the four-momenta of the produced charmed quarks. In our analysis below we 
will use the Sudakov decomposition, which has the following form: 

Pc = &lPe + PlPp + PcT, Pc = a 2 p e + P 2 Pp + PcT, ^ 

h = xtpe + k 1T , k 2 = x 2 p p + k 2T , 

where kyr, k 2 r, PcT and p e r are the transverse four-momenta of the corresponding particles. 
It is important that k^ T = — k\ T ^ and k| T = — k\ T ^ 0. If we make replacement k\ — > p e 
and set x\ = 1 and kyr = 0, then we easily obtain more simpler formulas corresponding to 
photon-gluon fusion subprocess. In the ep center-of-mass frame we can write 

p e = v^/2(l, 0,0,1), p p = v ^/2(l, 0,0,-1), (2) 

where s = (p e + p P ) 2 is the total energy of the process under consideration and we neglect 
the masses of the incoming particles. The Sudakov variables are expressed as follows: 

n ^cT ( v n m cT f v V 7 

Pi = exp(-y c ), j3 2 = —=r exp(-y g ), 

where m C T and m g T are the transverse masses of the produced quarks, and y c and y 5 are 
their rapidities (in the ep center-of-mass frame). From the conservation laws we can easily 
obtain the following conditions: 

x 1 = a l + a 2 , x 2 =[3 1 +[3 2 , k 1T + k 2T = p cT + p 5T . (4) 

In order to be sensitive to higher-order effects and to distinguish between direct-enriched and 
resolved-enriched regions the variable x° bs is often used [5] in the analysis of the data which 
contain the jets. This variable, which is the fraction of the photon momentum contributing 
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to the production of two jets with highest transverse energies E^, is experimentally defined 

as 

pjetl^-nJ 6 *! , 77ij ct 2 „-ni ot 2 

T obs _ E T e i +E T e 

where yE e is the initial photon energy and rf eti are the pseudo-rapidities of these hardest 
jets. The pseudo-rapidities rf cti are defined as rf eti = — lntan(^ eti /2), where 9 ieti are the 
polar angles of the jets with respect to the proton beam. The selection of x° bs > 0.75 and 
a;° bs < 0.75 yields samples enriched in direct and resolved photon processes, respectively. 



Another interesting variables, namely x° hs (D*) and cos9*, are also often used [5, 6] in the 
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analysis of the experimental data. The x° (D*) variable can be constructed in an analogous 
way to the traditional x° hs [5]. Using the D* ± meson and the untagged jet (i.e. jet which is 
not matched to a D*^ meson) of highest E^, the quantity x° hs (D*) is given by 

< s m = PTe ±|p , ( 5) 

where pt and rj are the transverse momentum and pseudo-rapidity of the produced D*^ 
meson. The scattering angle 9* is defined as 

cos0*=tanh- — . (11) 

Studying the distribution of cos 9* also gives us the possibility to learn about the size of the 
contributions from different production mechanisms [5]. Finally, the inelasticity z, defined 
by z = (p p ■ p) / (p p ■ q) with p and q being the four-momenta of the final D*^ meson and 
the exchanged photon, is a measure of the fraction of photon energy transferred to the D*^ 
meson in the proton rest frame. This quantity is sensitive to both the production mechanism 
and to the c — > D* fragmentation details [6] . 

2.2 Cross section for charm photoproduction 

The main formulas for the total and differential cross sections for charm production cross 
sections were obtained in our previous papers [11, 12, 26]. Here we recall some of them. In 
general case, the cross section a according to /^-factorization theorem can be written as a 
convolution 

a = J dk^a(k 2 T ,fi 2 )A(x,\4,fi 2 ), (6) 

where <5"(ky,/z 2 ) is the cross section corresponding to the relevant partonic subprocess under 
consideration and A(x, k^, fi 2 ) is the unintegrated gluon distribution. The direct photon 
contribution to the differential cross section of •yp — > cc + X process is given by 

da^^p^cc + X) f \M\ 2 ^g*^cc) 2 2 2 # 2 # c 

" dy c dp 2 cT ~ J 16tt(x 2 s) 2 (1 - ai) A W> k *r>l* ^ 2?r 2n , iV 

where \M\ 2 {ig* — > cc) is the squared off-shell matrix element which depends on the trans- 
verse momentum k 2T , 02 and C are the azimuthal angles of the initial virtual gluon and 
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the produced quark, respectively. The formula for the resolved photon contribution can be 
obtained by the similar way. But one should keep in mind that convolution in (6) should 
be made also with the unintegrated gluon distribution A^(x, k^,/i 2 ) in a photon. The final 
expression for the differential cross section has the form 

do-( res %p -> cc + X) _ f \M\ 2 {g*g* -> cc) 

■ X 



dy c dp 2 cT J 167r(xix 2 s) 2 , g , 

2 2 2 2 2 2 < ^'^ > " 

x Ay(xi, k 1T , n )A[x 2 , k^, // )<ik 1T <ik 2T <i?/e 



2tt 2tt 2tt 



where 0i is the azimuthal angle of the initial virtual gluon having fraction X\ of a initial 
photon longitudinal momentum. It is important that the squared off-shell matrix element 
\M\ 2 {g*g* — > cc) depends on the both transverse momenta k 2 T and k^. The analytic 
expressions for the |-M| 2 (7g* — > cc) and \M.\ 2 (g*g* — > cc) have been evaluated in our previous 
papers [12, 26]. Note that if we average (7) and (8) over kir and k 2 T and take the limit 
k 2 T — > and k?. T — > 0, then we obtain well-known formulas corresponding to the leading- 
order (LO) QCD calculations. 

The recent experimental data [5, 6] taken by the HI and ZEUS collaborations refer to 
the D*^ photoproduction in ep collisions, where the electron is scattered at small angle and 
the mediating photon is almost real (Q 2 ~ 0). Therefore the 'jp cross sections (7) and (8) 
need to be weighted with the photon flux in the electron: 

da(ep^cc + X) = J f 1 / e {y)dy da{^p — ► cc + X), (9) 

where y is a fraction of the initial electron energy taken by the photon in the laboratory 
frame, and we use the Weizacker- Williams approximation for the bremsstrahlung photon 
distribution from an electron: 

/<y/e ( y ) = ^ ( i + (i-y) 2 ln Cax , 2m 2 y M L\ V (io) 

Here a em is Sommerfeld's fine structure constant, m e is the electron mass, = m^y 2 / (1 — 
y) 2 and ~ 1 GeV 2 , which is a typical value for the recent photoproduction measure- 

ments at HERA. 

The multidimensional integration in (7), (8) and (9) has been performed by means of 
the Monte Carlo technique, using the routine Vegas [30]. The full C++ code is available 
from the authors on requeso This code is practically identical to that used in [11, 12], with 
exception that now we apply it to calculate inclusive and jet(s) associated charm production 
in another kinematical region. 



3 Numerical results 

We now are in a position to present our numerical results. First we describe our the- 
oretical input and the kinematical conditions. As it was mentioned above, the recent ex- 
perimental data [5, 6] on the charm photoproduction at HERA come from both HI and 

1 lipatov@theory.sinp. msu.ru 
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ZEUS collaboration. The ZEUS measurements [5] are performed in the following kinemat- 
ical region: 130 < W < 280 GeV, Q 2 < 1 GeV 2 , |^ ct | < 2.4, > 6 GeV, p T > 3 GeV 
and —1.5 < rj < 1.5. These data have been taken at a proton energy of 920 GeV and 
an electron energy of 27.5 GeV, which corresponds to a ep center-of-mass (cm.) energy of 
t/s = 318 GeV. Here and in the following all kinematic quantities are given in the laboratory 
frame where positive OZ axis direction is given by the proton beam. The more recent HI 
data [6] refer to the kinematical region defined by 171 < W < 256 GeV, Q 2 < 10~ 2 GeV 2 , 
|V et | < 1.5, E^ > 3 GeV, p T > 2 GeV and -1.5 < 77 < 1.5. 

3.1 Theoretical uncertainties 

There are several parameters which determined the overall normalization factor of the 
cross sections (7) and (8): the charm mass m c , the factorization and normalisation scales jip 
and ix r and the unintegrated gluon distributions in a proton A(x, kf,, fx 2 ) and in a photon 
v/4.ry(x, k^, fx ). 

Concerning the unintegrated gluon densities in a proton, we have tried three different 
sets of the unintegrated gluon densities in a proton, namely J2003 (set 1) [27], KMS [28] 
and KMR [29]. All these distributions are widely discussed in the literature (see, for exam- 
ple, review [31-33] for more information). Here we only shortly discuss their characteristic 
properties. First, the J2003 (set 1) gluon density has been obtained [27] from the numeri- 
cal solution of the full CCFM equation. The input parameters were fitted to describe the 
proton structure function F 2 (x,Q 2 ). Note that this density contain only singular terms in 
the CCFM splitting function P gg {z). The J2003 (set 1) distribution has been applied in 
the analysis of the forward jet production at HERA and charm and bottom production at 
Tevatron [27] (in the framework of Monte-Carlo generator Cascade [25]) and has been also 
used in our calculations [11, 12, 26]. 

Another set (the KMS) [28] was obtained from a unified BFKL-DGLAP description of 
early F 2 (x ) Q 2 ) data and includes the so-called consistency constraint [34]. The consistency 
constraint introduces a large correction to the LO BFKL equation. It was argued [34] that 
about 70% of the full NLO corrections to the BFKL exponent A are effectively included in 
this constraint. The KMS gluon density is successful in description of the beauty hadropro- 
duction at Tevatron [15] and photoproduction at HERA [12]. 

The last, third unintegrated gluon distribution A(x, k^, fx 2 ) used here (the so-called KMR 
distribution) is the one which was originally proposed in [29]. The KMR approach is the 
formalism to construct unintegrated gluon distribution from the known conventional parton 
(quark and gluon) densities. It accounts for the angular-ordering (which comes from the 
coherence effects in gluon emission) as well as the main part of the collinear higher-order 
QCD corrections. The key observation here is that the fx dependence of the unintegrated 
parton distribution enters at the last step of the evolution, and therefore single scale evolution 
equations (DGLAP or unified BFKL-DGLAP) can be used up to this step. Also it was 
shown [29] that the unintegrated distributions obtained via unified BFKL-DGLAP evolution 
are rather similar to those based on the pure DGLAP equations. It is because the condition 
of the angular ordering constraint is more important [29] than including the BFKL effects. 
Based on this point, in the present paper we use much more simpler DGLAP equation up 
to the last evolution step in the case of application of the KMR procedure. In the numerical 
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calculations we have used the standard GRV (LO) parametrizations [35] of the collinear quark 
and gluon densities. Note that the KMR unintegrated parton distributions in a proton were 
used, in particular, to describe the prompt photon photoproduction at HERA [13, 36] and 
prompt photon hadroproduction Tevatron [14, 37]. 

In the case of a real photon, we have tested two different sets of the unintegrated gluon 
densities A^(x, k 2 ^, /i 2 ). First of them was obtained [38] from the numerical solution of the 
full CCFM equation (which has been also formulated for the photon). Here we will use 
this gluon density together with the J2003 (set 1) distribution when calculating the resolved 
photon contribution (8). Also in order to obtain the unintegrated gluon density in a photon 
we will apply the KMR method [29] to the standard LO GRV parton distributions [35]. In 
the numerical calculations we will use it together with the KMR distributions in a proton. 
Note that both gluon densities »4 7 ( ) discussed here have been already applied in 

the analysis of the charm and beauty quark [26, 39] and J meson production [26] in 77 
collisions at LEP2. 

We would like to point out that at present there is no the unintegrated gluon distribution 
corresponding to the unified BFKL-DGLAP evolution in a photon. Therefore the resolved 
photon contribution (8) is not taken into account in the case of KMS gluon distribution. 

Significant theoretical uncertainties in our results are connected with the choice of the 
factorization and renormalization scales. The first of them is related to the evolution of the 
gluon distributions, the other is responsible for the strong coupling constant cc s (//fj). As 
it often done for charm production, we choose the renormalization and factorization scales 
to be equal: Hr = fip = jj, = i\Jm1 + (p 2 -), where (p 2 -.) is set to the average p 2 - of the 
charmed quark and antiquark. In order to investigate the scale dependence of our results 
we will vary the scale parameter £ between 1/2 and 2 about the default value £ = 1. Note 
that we use special choice /i 2 = in the case of KMS gluon, as it was originally proposed 
in [27]. The fragmentation c — > D* is described by Peterson fragmentation function [40] 
with e c = 0.035 [41]. The branching ratio f(c — > D*) was set to the value measured by the 
OPAL collaboration: f(c — > D*) = 0.235 [42]. For completeness, we take the charmed quark 
mass m c = 1.4 GeV and use LO formula for the coupling constant a s (/i 2 ) with rif = 4 active 
quark flavours at A QCD = 200 MeV, such that a s (M 2 z ) = 0.1232. 

3.2 Inclusive D* ± production 

The results of our calculations are shown in Figs. 1 — 4 in comparison to the recent HI 
experimental data [6]. Instead of presenting our theoretical predictions as continuous lines, 
we adopt here the binning pattern encoded in the experimental data. The solid, dashed and 
dash-dotted histograms correspond to the results obtained with the J2003 (set 1), KMR and 
KMS unintegrated gluon densities, respectively. We observe a reasonable good agreement 
between our predictions and the HI experimental data [6] in most of the bins. However, there 
are some insignificant discrepancies. So, the calculated transverse momentum distribution 
falls less steeply at large pt than it is visible in the data. The similar effect was pointed 
out also in [6] where the massive NLO pQCD approach [18] (fmnr program [43]) and the 
Monte-Carlo generator cascade [35] have been used. Concerning the dependence of our 
predictions on the evolution scheme, we have found that the D*^ transverse momentum 
distribution is only little sensitive to the choice of unintegrated gluon density: very similar 
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predictions are obtained by all parametrizations (see Fig. 1). This is in agreement with 
results of the previous investigations [23]. The same situation is observed in W distribution 
where all gluon densities under consideration predict a very similar shapes. In contrast, the 
calculated pseudo-rapidity and z distributions strongly depend on the unintegrated gluon 
density used. Note that pseudo-rapidity distribution additionally has been measured [6] 
in the three bins of pt (see Fig. 2). In all these kinematical regions KMR gluon tends to 
overestimate the data in the rear direction and KMS tends to underestimate the data in the 
forward one. This is, in particular, due to the fact that gluon-gluon fusion (resolved photon 
contribution) is not taken into account in the case of KMS gluon in our calculations. At 
the same time the pseudo-rapidity distributions obtained with the J2003 (set 1) and KMR 
unintegrated gluon densities (where gluon-gluon fusion contribution is included everywhere) 
agree well with the HI data [6] in the forward pseudo-rapidity region. This is in agreement 
with the general expectations that resolved photon contributions are important at r\ > 0. 

We would like to point out also that gluon-gluon fusion contributes significantly in the 
low z region. In order to illustrate this effect we have separately shown the contributions 
from the photon-gluon (dashed histogram) and gluon-gluon fusion (dash-dotted histogram) 
mechanisms within the kinematic range of the HI measurement [6] (see Fig. 3). The solid 
histogram represents the sum of both these contributions. We have used here the KMR 
unintegrated gluon density. It is clear that gluon-gluon fusion mechanism is important at 
low z and should be taken into account in description of the experimental data. 

Next, the total inclusive -D* ± photoproduction cross section a(ep — > e'D*^ + X) has 
been measured [6] and it was found to be equal to 6.45 ± 0.46 (stat.) ±0.69 (sys.) nb. The 
results of our calculations supplemented with the different unintegrated gluon densities are 
collected in Table 1. The theoretical uncertainties (which are given for the J2003 (set 1) 
and KMR distributions) are connected with variation on the scale /i 2 as it was described 
above. The predictions of cascade, pythia as well as NLO pQCD calculations (GMVFNS 
approach and FMNR program) are shown for comparison. The central values from massive 
NLO pQCD calculations (fmnr program) and from cascade are slightly lower than the 
measured result, whereas those of pythia [44] and GMVFNS [20] are higher. One can see 
that our predictions are rather close to ones from the cascade and fmnr programs and 
reasonably agree with the HI data within the theoretical and experimental uncertainties. 
Also we estimate the individual contributions from the photon-gluon and gluon-gluon fusion 
to the total cross section in the /^-factorization approach. We have found it to be about 80% 
and 20%, respectively. Additionally we investigate the scale dependence of our predictions 
(see Fig. 4). In these plots the solid histograms were obtained by fixing both the factorization 
and normalization scales at the default value fi 2 , whereas upper and lower dashed histograms 
correspond to the scale variation as it was described above. Here we have used the KMR 
unintegrated gluon density. One can see that scale variation changes the normalization of 
predicted cross section by 20 - 30% approximately. Note that we have not varied the charmed 
quark mass and used the default value of m c in all calculations. 

3.3 Associated D* ± and single jet production 

Now we demonstrate how /c T -factorization approach can be used to calculate the semi- 
inclusive charm photoproduction rates. The basic photon-gluon or gluon-gluon fusion sub- 
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Source 


a(ep -> e'D** + X) [nb] 


HI measurement [6] 

CASCADE [25] 
PYTHIA [45] 

rMINK. r±^\ 

GMVFNS [20] 
J2003 (set 1) 
KMR 
KMS 


6.45 ± 0.46 (stat.) ± 0.69 (sys.) 

r OO+0.54 
o -< 3 °-0.62 

8.9 

q o+5.3 
0-^-4.0 

4.921^ 
6.57±;;S 
5.10 



Table 1: The total cross section of the inclusive D*^ photoproduction in electron-proton 
collisions at Q 2 < 10~ 2 GeV 2 , 0.29 < y < 0.65, p T > 2 GeV and -1.5 < rj < 1.5. 

processes under consideration give rise to two high-energy charmed quarks, which can further 
evolve into hadron jets. In our calculations the produced quarks (with their known kine- 
matical parameters) were taken to play the role of the final jets. These two quarks are 
accompanied by a number of gluons radiated in the course of the gluon evolution. As it 
has been noted in [22], on the average the gluon transverse momentum decreases from the 
hard interaction block towards the proton. As an approximation, we assume that the gluon 
emitted in the last evolution step and having the four-momenta k' compensates the whole 
transverse momentum of the gluon participating in the hard subprocess, i.e. k' T = —k T . All 
the other emitted gluons are collected together in the proton remnant, which is assumed to 
carry only a negligible transverse momentum compared to k'^. This gluon gives rise to a final 
hadron jet with = [k'^l in addition to the jet produced in the hard subprocess. From 
these three hadron jets we choose the one jet (or two jets) carrying the largest transverse 
energy, and then compute the charm and associated jet(s) production rates. 

In the recent analysis [5, 6] performed by of the HI and ZEUS collaborations jets are 
divided into two categories: jets of the first category are associated with the D*^ meson (D*- 
tagged jet), while jets of the second category are not matched to a D*^ meson (/}*-untagged 
jet). The inclusive, _D*-tagged and untagged jet cross sections have been measured [5] by the 
ZEUS collaboration whereas HI data refer to the D* and untagged jet cross sections only [6]. 
Note that formulation of the ZMVFNS do not provide the possibility to simultaneously 
determine kinematic variables related to the meson and the jet that includes it, so that 
a comparison to the massless NLO pQCD calculations is only available for the untagged jet 
cross sections [45]. In the following we will systematically compare the predictions from the 
^-factorization approach to all published data on associated charm and single jet production 
at HERA. 

3.3.1 D*-tagged jets: pj? and rf ct distributions 

The results of our calculations are shown in Fig. 5 in comparison to the ZEUS experimen- 
tal data [5]. Notations of all histograms here are the same as in Fig. 1. The pseudo-rapidity 
if et distribution additionally has been measured in the three bins of E^. One can see that 
shape of the data is well described by our predictions, although the normalisation is un- 
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derestimated (by a factor of 1.5). However, this discrepancy is not dramatic, because some 
reasonable variations in charm mass m c , energy scale fi 2 or Aqcd parameter (not shown in 
Figs.) can partially cover the visible disagreement. It is interesting that difference between 
the theoretical predictions calculated with different unintegrated gluon densities in a proton 
are somewhat less pronounced in comparison to the inclusive D*^ production case. The 
KMS gluon tends to predict a larger cross sections than ones obtained with other gluon 
densities under consideration. This fact is in contrast with the results presented in the pre- 
vious section where KMR gluon dominates. The possible explanation is connected with the 
difference in the kinematical region (in which our predictions as well as experimental data 
were presented). We would like to point out also that our central predictions are very similar 
to the massive NLO pQCD results [18]. It was claimed [5] that normalisation of the ZEUS 
data for all distributions is reasonable described by the upper limit of these NLO pQCD 
calculations. 

3.3.2 /^-untagged jets: x° hs (D*) distribution 

Now we turn to the D* ± and untagged jet production. First we discuss the very inter- 
esting subject of study which is connected with the individual contributions from the direct 
and resolved photon mechanisms. As it was already mentioned above, the £° bs (.D*) variable 
(which corresponds at leading order to the fraction of the exchanged photon momentum in 
the hard scattering process) provides a tool to investigate the relative importance of these 
different contributions. In LO collinear approximation, direct photon events at parton level 
have x^ s (D*) ~ 1, while the resolved photon events populate the low values of a;° bs (£>*). 
The same situation is observed in a NLO calculations, because in the three parton final 
state any of these partons are allowed to take any kinematically accessible value. In the 
/cT-factorization formalism the hardest transverse momentum parton emission can be any- 
where in the evolution chain, and does not need to be closest to the photon as required 
by the strong /x 2 ordering in DGLAP. Thus, if hardest jet originates from the cc pair, then 
a;° bs (-D*) is close to unity, but if a gluon from the initial cascade form the hardest trans- 
verse momentum jet, then x^ s (D*) < 1. This statement is clearly demonstrated in Fig. 6 
where separately shown the contributions from the photon-gluon (dashed histogram) and 
gluon-gluon fusion (dash-dotted histogram) subprocesses within the kinematic range of the 
ZEUS measurement [5]. The solid histogram represents the sum of both these contributions. 
We have used here the KMR unintegrated gluon density for illustration. In agreement with 
the expectation for direct photon processes, the peak at high values of the x° bs (.D*) is ob- 
served. However, one can see that off-shell photon-gluon fusion results also in substantial 
tail at small values of x° bs (D*). The existence of this plateau in the collinear approximation 
of QCD usually is attributted to the charmed quark excitation from resolved photon. In 
the /cy-factorization approach such plateau indicates the fact that the gluon radiated from 
evolution cascade appears to be harder than charmed quarks (produced in hard parton in- 
teraction) in a significant fraction of events [11, 22]. Therefore we can conclude that the 
/c^-factorization formalism effectively imitates the charm component of the photon [11, 22]. 
However, the predicted tail at small x° hs (D*) values is strongly depends on the unintegrated 
gluon distributions used, as it was demonstrated in [11]. Note that the gluon-gluon fusion 
events (with a gluon coming from the photon) are distributed over the whole £° bs (-D*) range. 
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It is clear that these events play role at small values of a;° bs (.D*) only and that now the total 
contribution to the cross section from the gluon inside the resolved photon is too small to 
be useful^. This fact is in contrast with the inclusive D*^ production case and coincide with 
the results [45] obtained in the massless NLO QCD approximation. 

In Fig. 7 we confront the calculated x° bs (.D*) distribution with the ZEUS data [5]. One 
can see that our results corresponding to different gluon densities do not agree well with 
the ZEUS data. If this disagreement is true, then the possible explanation of this fact is 
the following. The calculated cross sections at low x^ s (D*) are defined by the average 
value of the gluon transverse momenta (hx) which is generated in the course of the non- 
colliner evolution. It is due to the fact that events when the gluon jet has the largest pt 
among the three hadron jets contribute only in this kinematical region. So, this average 
gluon (h?) which generated by the all three versions of the unintegrated gluon distributions 
under consideration possible could be too small to describe the ZEUS data. The similar 
situation was observed in [11] where the D*^ and dijet associated photoproduction has been 
considered. However, the /cy-factorization approach well describes the experimental data 
with the cut on the dijet invariant mass M > 18 GeV [11]. It demonstrates that this 
cut is essential for applicability of the description of resolved photon contributions by non- 
collinear evolution in a proton only. Therefore we can conclude that further theoretical and 
experimental investigations are necessary in order to understand and adequatively describe 
the ZEUS experimental data in the low x^ s (D*) region. 

3.3.3 D*-untagged jets: p T , rj, pj? and rf et distributions 

In Figs. 8 — 10 we confront the differential cross section of ep — ► D* ± + jet + X in 
photoproduction as measured [5, 6] by HI and ZEUS collaborations with our predictions. 
One can see in Fig. 8 that the differential cross section as a function of r] and p^ measured 
by the HI collaboration is well described (both in shape and in magnitude) by our calcula- 
tions. Also the predicted cross section as a function of r\ — r/J et (see Fig. 10) are in a good 
agreement with the experimental data. Note that there is some difference in shape of the 
r]i et distribution between the data and our calculations. Similar to inclusive D*^ production 
case, all unintegrated gluon densities under consideration slightly overestimate the da/dpx 
distribution at large p? (see Fig. 8). But in general the HI data on the total and differential 
cross sections are reasonable well reproduced by our calculations. In Table 2 the estimations 
of total cross section obtained from our calculations as well as from cascade, pythia, 
massive and massless NLO pQCD evaluations are listed and compared with the HI data [6]. 
The theoretical uncertainties (which are given for the J2003 (set 1) and KMR distributions) 
are connected with variation on the scale /i 2 , as it was made earlier in the inclusive D*^ 
production case. 

Concerning the ZEUS measurement, the situation is not clear, as it seen from Fig. 9. 
One can see that transverse energy distribution is underestimated in all bins. The 
difference in normalisation between our predictions and the ZEUS data is about 1.5 or 

2 This is one of reasons why the values of our cross sections (see Table 1 and 2) obtained with the 
J2003 (set 1) unintegrated gluon distribution are not more the Cascade predictions. The main reasons of 
this distinction are the different factorization scales in our calculations and the Cascade ones and also the 
account for the initial and final showers in the MC generator Cascade [6]. 
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Source 


a(ep -> e'D** + jet + X) [nb] 


HI measurement [6] 

CASCADE [25] 
PYTHIA [45] 

rlVlINrt n"| 

ZMVFNS [46] 
J2003 (set 1) 
KMR 
KMS 


3.01 ± 0.29 (stat.) ± 0.33 (sys.) 
3.081°;! 
3.8 

z - oo -0.42 
o np; +0.62 

o a 7+O.6O 
-0.39 
o fi 1+0.84 
°- ui -0.68 

3.11 



Table 2: The total cross section of the associated D* ± and single jet photoproduction in 
electron-proton collisions at Q 2 < 0.01 GeV 2 , 0.29 < y < 0.65, p T > 2 GeV, -1.5 < rj < 1.5, 
p£* > 3 GeV and -1.5 < ?? jct < 1.5. 

even 2. This fact is in clear contrast with comparison of our theoretical results with the 
HI measurements. Also the shape of measured rf et distributions is very different from the 
calculated one. Particularly, the ZEUS data are overshoot our theoretical estimations in 
the forward direction. But a good agreement is achieved in the rear pseudo-rapidity region 
(namely, at rf ct < 0.5). The similar situation is obtained [5] in the framework of both massive 
and massless NLO pQCD approaches [43, 45]. 

Now we turn to some qualitative comparison between the predictions of /^-factorization 
approach and the ones obtained in the framework of collinear NLO pQCD approximation. 
It is well known that in the case of inclusive single hadron photoproduction the direct 
and resolved photon contributions are accumulated in the backward and forward directions, 
respectively. Using the ZMVFNS scheme, it was shown [45] that the pseudo-rapidity rf et of 
.D*-untagged jet (but not pseudo-rapidity 77 of the meson) can serve as a discriminator 
between these two contributions. It is because the shapes of direct and resolved photon 
mechanism are very different to each other. In Fig. 11 we split our central predictions for 
da/drf ct distribution calculated in the ZEUS kinematical region into their direct and resolved 
photon components using the a;° bs (-D*) variable and show the results as the solid and dashed 
histograms, respectively. We have used the KMR unintegrated gluon densities in a proton 
and in a photon. Similar to massless NLO pQCD [45] , we obtain a strong different behaviour 
of calculated cross sections at x^ s (D*) > 0.75 and x^ s (D*) < 0.75. Moreover, we reproduce 
well the direct contribution both in normalization and shape (see also Fig. 8 from Ref. [45]). 
But predicted resolved photon component lie below NLO pQCD one by a factor of about 2, 
although agree with NLO pQCD results in a shape. The observed difference in normalisation 
is strongly depends, of course, on the unintegrated gluon distribution used. In particular, 
difference between our calculations and results presented in [45] indicates again that average 
gluon (kr) is too small and that further theoretical investigations of non-collinear gluon 
evolution are needed. However, at qualitative level, the //^-factorization formalism in a 
simplest way reproduces well the basic properties of the standard fixed-order calculations 
for process under consideration. 
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3.3.4 D*-untagged jets: cos 9* variable 

Next we turn to the different angular correlations in the charm production at HERA. First 
we discuss the cos 9* variable. As it was already mentioned above, studying of distribution 
on cos 9* also give us the possibility to learn about the size of the contribution from different 
production mechanisms. It is because this quantity is sensitive to the spin of the propagator 
in the hard subprocess [5]. In direct photon processes — > cc the propagator in the 
LO QCD diagrams is a quark, and the differential cross section rises slowly towards high 
| cos 6 1 * | values, namely proportional to (1 — |cos^*|) _1 . In resolved processes, the gluon 
propagator is allowed at LO and dominates over the quark propagator due to the stronger 
gluon-gluon coupling compared to the quark-gluon coupling. If most of the resolved photon 
events are produced as a result of charm from the photon, a gluon-exchange contribution 
in eg — > eg subprocess should dominate. This results in a steep rise of the cross section 
towards high |cos#*| values: da/dcos9* ~ (1 — |cos#*|) -2 . A recent ZEUS analysis [5] on 
dijet angular distributions in the D* ± photoproduction has shown that the measured cross 
section from resolved-enriched events (i.e. events with x^ s (D*) < 0.75) exhibits a distinct 
asymmetry with a strong rise towards cos9* = —1. This behaviour suggest that events with 
£° bs (.D*) < 0.75 are dominantly produced by charmed quarks coming from the photon side. 
On the other hand, the cos#* distribution for direct-enriched events (where x^ s (D*) > 0.75) 
is almost symmetric. In the /c^-factorization approach the cos#* distribution is determined 
only by the photon-gluon fusion off-shell matrix element which cover both scattering process 
(since there is no restriction on the transverse momenta along the evolution cascade, as it was 
already discussed above). In order to study this quantity in a detail, in Fig. 12 we show the 
differential cross section as a function of cos 9* for the direct-enriched (solid histogram) and 
resolved-enriched (dashed histogram) samples separately within the kinematic range of the 
ZEUS experiment. We have used here the KMR unintegrated gluon density for illustration. 
One can see that resolved photon-like events exhibit a strong rise towards cos#* = —1, i.e. 
in photon direction, consistent with a dominant contribution from gluon exchange. In our 
theoretical calculations, the peak at cos#* = —1 at low x^ s (D*) clearly illustrates again 
that the /^-factorization approach effectively reproduces the charm excitation processes 
using only the photon-gluon fusion off-mass shell matrix elements. However, we should 
point out that the absolute normalization of this peak differs from the one calculated [45] 
in the massless NLO pQCD approximation (by a factor of about 2). This fact is in a 
full agreement with the results presented in the previous section (in Fig. 11) and indicates 
again that resolved photon contribution is underestimated by our calculations (compared to 
the NLO pQCD ones). Unfortunately, at present there is no experimental data on cos 9* 
distribution in the associated and single jet photoproduction. Therefore additional 
experimental efforts in this field possible can give us the possibility to better constraint the 
unintegrated gluon density in a proton. 

3.3.5 D*-untagged jets: azimuthal correlations 

Further understanding of the process dynamics and in particular of the high-order cor- 
rection effects may be obtained from the transverse correlation between the produced D*^ 
meson and the jet. The HI collaboration has measured [6] the distribution on the A0 angle, 
which is the difference in azimuth between the D*^ meson and the jet. In the naive leading 
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order approximation, the distribution over A<j) must be simply a delta function 5(A<f) — ir), 
since the produced Z)* ± meson and the jet are back-to-back in the transverse plane. Taking 
into account the non-vanishing initial parton transverse momenta k 1T and k 2 r leads to the 
violation of this back-to-back kinematics in the /c^-factorization approach [15]. 

The calculated A<p distributions are shown in Fig. 13 in comparison to the HI data. Addi- 
tionally, in Figs. 14 and 15 we plot our predictions separately for the regions x^ s (D*) > 0.75 
and a;° bs (.D*) < 0.75, where direct and resolved photon induced processes dominate, respec- 
tively. Last results have been calculated for the ZEUS kinematical range, as it was described 
above. It is clear that /c T -factorization approach predicts a large fraction of events where the 
D*^ and the jet are not back-to-back. The measured deviations from the back-to-back topol- 
ogy are reasonable well described by our calculations. Best agreement is achieved using the 
KMR unintegrated gluon density. There is general behaviour of theoretical results obtained 
using the unintegrated gluon densities under consideration. All of them overestimate the HI 
data at A0 ~ (where gluon k T is large) as well as at A0 ~ n, i.e. in the back-to-back 
region. But at moderate A0 our predictions tends to underestimate the data. One can see 
that region of low A<p is very sensitive to the unintegrated gluon distributions: the difference 
between predictions can be one order of magnitude (see Figs. 14 and 15). The similar results 
has been obtained in [11, 15]. Therefore we can expect that further theoretical and exper- 
imental studying of these correlations will give important information about non-collinear 
parton evolution dynamics in a proton and in a photon. 

4 Conclusions 

We have investigated the charm photoproduction in electron-proton collisions at HERA 
using the /cy-factorization approach of QCD. Our analysis cover the inclusive charm produc- 
tion as well as charm and associated jet production processes. Both photon-gluon (direct) 
and gluon-gluon fusion (resolved) mechanisms were taken into account. Our investigations 
were based on the leading-order off-shell matrix elements of the relevant partonic subpro- 
cesses. The total and differential cross sections have been calculated and the comparisons to 
the recent HI and ZEUS experimental data have been made. In numerical analysis we have 
used the unintegrated gluon densities which are obtained from the full CCFM, from unified 
BFKL-DGLAP evolution equations as well as from the Kimber-Martin-Ryskin prescription. 

We have shown that the /c^-factorization approach reproduces reasonably well the nu- 
merous HI data on both inclusive charm and associated charm and jet production. At the 
same time we find that the ZEUS measurements overshoot our theoretical predictions by a 
factor of about 1.5 or even 2. It means that additional efforts should be done from both the 
theoretical and experimental sides in order to reduce this disagreement. 

We find that our results agree with the standard NLO pQCD ones at qualitative level 
for process under consideration. It was demonstrated that off-shell matrix elements com- 
bined with the non-collinear evolution of gluon densities in a proton effectively simulate 
the charmed quark excitation processes. We find that resolved photon contributions are 
underestimated in the framework of /c T -factorization formalism compared to the massless 
NLO pQCD calculations. Of course, degree of this underestimation strongly depends on the 
unintegrated gluon density used. 
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Special attention in our calculations has been drawn to the specific angular correlations 
between the produced D*^ 1 meson and jets in final state. In particular, we demostrate the 
strong sensitivity of the A0 distribution at low A<p to the unintegrated gluon densities in a 
proton and in a photon. 

To conclude, we believe that the /^-factorization formalism holds a possible key to un- 
derstanding charm production at HERA. However, there are still large uncertainties, and 
much more work needs to be done. 
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Figure 1: The differential cross sections of the inclusive D*^ production at HERA calculated 
in the kinematic range Q 2 < 0.01 GeV 2 , 0.29 < y < 0.65, p T > 2 GeV and -1.5 < i] < 
1.5. The solid, dashed and dash-dotted histograms correspond to the predictions obtained 
with the J2003 (set 1), KMR and KMS unintegrated gluon densities, respectively. The 
experimental data are from HI [6]. 
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Figure 2: Inclusive D* ± cross sections as a function of rj for three bins of p T - Notations of 
all histograms here are the same as in Fig. 1. The experimental data are from HI [6]. 
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Figure 3: The da/dz distribution of the inclusive D*^ production at HERA calculated in 
the kinematic range Q 2 < 0.01 GeV 2 , 0.29 < y < 0.65, p T > 2 GeV and —1.5 < rj < 
1.5. Separately shown the contributions from the photon-gluon (dashed histogram) and 
gluon-gluon (dash-dotted histogram). Solid histogram represents the sum of both these 
contributions. The KMR unintegrated gluon densities in a proton and in a photon have 
been used. 
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Figure 4: The da/dpT distribution of the inclusive D* ± production at HERA calculated in 
the kinematic range Q 2 < 0.01 GeV 2 , 0.29 < y < 0.65 and -1.5 < r) < 1.5. Solid, upper 
dashed and lower dashed histograms were obtained with the /i 2 = rri^, fj, 2 = l/Am? and 
ji 2 = 4m 2 -, respectively. The experimental data are from HI [6]. 
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Figure 5: The differential cross sections of the .D and tagged jet production at HERA 
calculated in the kinematic range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, 
\r)\ < 1.5, —1.5 < rf ct < 2.4 and > 6 GeV. Notations of all histograms here are the same 
as in Fig. 1. The experimental data are from ZEUS [5]. 
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Figure 6: Direct (dashed histogram) and resolved photon (dashed-dotted histogram) com- 
ponents of our central predictions of the differential cross section da / dx^ s (D*) calculated 
in the kinematic range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, \r]\ < 1.5, 
— 1.5 < r] sct < 2.4 and E 3 ^ > 6 GeV. The KMR unintegrated gluon densities in a proton and 
in a photon have been used. 
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Figure 7: The da / dx'^ >s {D*) distribution of the associated D*^ and jet production at HERA 
calculated in the kinematic range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, \rj\ < 1.5, 
— 1.5 < 77 jct < 2.4 and > 6 GeV. Notations of all histograms here are the same as in 
Fig. 1. The experimental data are from ZEUS [5]. 
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Figure 8: The differential cross sections of the and untagged jet production at HERA 
calculated in the kinematic range Q 2 < 0.01 GeV 2 , 0.29 < y < 0.65, p T > 2 GeV, -1.5 < 
rj < 1.5, J?? > 3 GeV and —1.5 < rf ct < 1.5. Notations of all histograms here are the same 
as in Fig. 1. The experimental data are from HI [6]. 
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Figure 9: The differential cross sections of the and untagged jet production at HERA 
calculated in the kinematic range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, \r}\ < 1.5, 
— 1.5 < rf ct < 2.4 and > 6 GeV. Notations of all histograms here are the same as in 
Fig. I. The experimental data are from ZEUS [5]. 
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Figure 10: The cross section of the D*^ and untagged jet production as a function of 
i] — rf ct calculated in the kinematic range Q 2 < 0.01 GeV 2 , 0.29 < y < 0.65, p T > 2 GeV, 
—1.5 < r] < 1.5, > 3 GeV and —1.5 < rf ct < 1.5. Notations of all histograms here are 
the same as in Fig. 1. The experimental data are from HI [6]. 



28 



1.0 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



0.8 - 



9 0.6 - 



0.4 



0.2 - 







i — r— 'i | | | | | | | | | | | | | | | l 



-1 



1 

n iet 



Figure 11: Direct (solid histogram) and resolved photon (dashed histogram) components of 
our central predictions of the differential cross section da/drj ict calculated in the kinematic 
range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, \r]\ < 1.5, -1.5 < ?? jct < 2.4 and 
> 6 GeV. The KMR unintegrated gluon densities in a proton and in a photon have been 
used. 
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Figure 12: Direct (solid histogram) and resolved photon (dashed histogram) components of 
our central predictions of the differential cross section da/d cos 9* calculated in the kinematic 
range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, \t]\ < 1.5, -1.5 < ?? jct < 2.4 and 
> 6 GeV. The KMR unintegrated gluon densities in a proton and in a photon have been 
used. 
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Figure 13: The cross section of the D*^ and untagged jet production as a function of 
A0 calculated in the kinematic range Q 2 < 0.01 GeV 2 , 0.29 < y < 0.65, p T > 2 GeV, 
— 1.5 < rj < 1.5, pj? > 3 GeV and —1.5 < rf ct < 1.5. Notations of all histograms here are 
the same as in Fig. 1. The experimental data are from HI [6]. 
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Figure 14: The cross section of the D*^ and untagged jet production as a function of A(f> 
calculated in the kinematic range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, \rj\ < 1.5, 
-1.5 < ?7 jet < 2.4, > 6 GeV and x° hs (D*) > 0.75. Notations of all histograms here are 
the same as in Fig. 1. 
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Figure 15: The cross section of the D*^ and untagged jet production as a function of A(f> 
calculated in the kinematic range Q 2 < 1 GeV 2 , 130 < W < 280 GeV, p T > 3 GeV, \rj\ < 1.5, 
-1.5 < ?7 jet < 2.4, > 6 GeV and x° hs (D*) < 0.75. Notations of all histograms here are 
the same as in Fig. 1. 
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